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We analyse a problem of anti -plane shear in a bi-material plane containing a semi-infinite crack situated on a soft imperfect interface. 
The plane also contains a small thin inclusion (for instance an ellipse with high eccentricity) whose influence on the propagation of 
the main crack we investigate. An important element of our approach is the derivation of a new weight function (a special solution to 
a homogeneous boundary value problem) in the imperfect interface setting. The weight function is derived using Fourier transform 
and Wiener-Hopf techniques and allows us to obtain an expression for an important constant cr (which may be used in a fracture 
criterion) that describes the leading order of tractions near the crack tip for the unperturbed problem. We present computations that 
demonstrate how <t varies depending on the extent of interface imperfection and contrast in material stiffness. We then perform 
perturbation analysis to derive an expression for the change in the leading order of tractions near the tip of the main crack induced 
by the presence of the small defect, whose sign can be interpreted as the inclusion's presence having an amplifying or shielding 
effect on the propagation of the main crack. 
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conditions, whose derivation can be found for example in An 



1. Introduction 

In this paper we present a method to evaluate important con- 
stants which describe the behaviour of physical fields near crack 
tips in a perturbed problem set in a domain containing an im- 
perfect interface. 

Imperfect interfaces account for the fact that the interface be- 
tween two materials is almost never sharp. Atk inson| ( |1977| l ac- 
counted for this observation by by placing a very thin strip of a 
homogeneous material in the model between two larger bodies 
with different elastic moduli to that of the strip. If the thin layer 
is considered to be either much softer or stiffer than the main 
bodies, its presence can be replaced in models by transmission 



X 



tipov et al. (2001 1 for a soft imperfect interface, or Mishuris et 
al. (2006) for a stiff imperfect interface. We shall consider only 
soft imperfect interfaces in the present paper. 



Klarbring et al. ( 1998 1 presented an asymptotic model of ad- 



hesive joints in a layered structure. Mishuris (2001 1 found the 
asymptotic behaviour of displacements and stresses in a vicinity 
of the crack tip situated on a soft imperfect interface between 
two different elastic materials, where the non-ideal interface is 
replaced by non-ideal transmission conditions. For such a case, 
the asymptotics are of a markedly different form to the perfect 
interface case, in which components of stress exhibit a square 
root singularity at the crack tip; such behaviour is not present 
for imperfect interface cracks. 

A key element of our approach will be the derivation of a 
new weight function. The concept of weight functions was in- 



troduced by |Bueckner| ( |1970] l. In the perfect interface setting 
these provide weights for the loads applied to the crack surfaces 
such that their weighted integrals over the crack surfaces pro- 



vide the stress intensity factors at a certain point. Vellender et 
[aT1 ( |201 1) modified the weight function technique to yield sim- 
ilarly useful asymptotic constants that characterise stress fields 
near crack tips along an imperfect interface. 

A survey of macro-microcrack interaction problems can be 
found in |Petrova et al.| ( 2000 1. Of particular relevance is the 
recent manuscript of Mishur is et al.| ( |20"TT| which examines an 
analagous problem to that presently considered with a perfect 
interface in place of the imperfect interface. The approach in 



that paper utilises the dipole matrix approach of Movchan et al. 
( |1995| l to construct an asymptotic solution that takes into ac- 
count the presence of a micro-defect such as a small inclusion. 
The present paper seeks to adapt this approach to the imperfect 
interface setting. 

2. Structure and summary of main results 

We adopt the following structure for the paper. We first for- 
mulate the physical problem before giving the weight function 
problem formulation. Fourier transform techniques allow us 
to obtain a Wiener-Hopf type problem for the weight func- 
tion, whose kernel we factorise in a computationally convenient 
fashion. The Wiener-Hopf equation is solved to yield expres- 
sions for the weight function and comparisons are drawn be- 
tween the perfect and imperfect interface weight function prob- 
lems. 
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We then use the reciprocal theorem (Betti formula) in the 
spirit of Willis et al. ( 1995] > to relate the sought physical solu- 
tion to the weight function. The presence of imperfect interface 
transmission conditions alters properties of the functions in the 
Betti identity and so different analysis is required. The applica- 
tion of Betti's identity enables us to find an expression for the 
leading order of tractions <x near the crack tip in terms of the 
new weight function and the imposed arbitrary tractions pre- 
scribed on the faces of the crack: 
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Here, bars denote Fourier transform, is a constant depend- 
ing on the material parameters and extent of interface imper- 
fection, [[ I/] and (IT) are respectively the jump and average of 
the weight function across the crack/interface line, and [p] and 
(p) are the jump and average of the tractions prescribed on the 
crack faces. 

In Section [7] we perform perturbation analysis to determine 
the impact on the tractions near the crack tip of the presence of 
a small inclusion. The asymptotic solution is sought in the form 



u(x, e) = u m {x) + £W u, (f) + e l u w {x) + o(e l ), e -> 0, (2) 

where m (0) is the unperturbed physical displacement solution 
(the solution with no inclusion present), sW is a boundary layer 
concentrated near the inclusion and s 2 m (1) is introduced to fulfil 
the original boundary conditions on the crack faces and along 
the imperfect interface. This enables us to find the first order 
variation in the crack tip tractions; we expand the constant ctq 
as 

o-o = erf + e 2 Acr + o(e 2 ), s -> 0, (3) 
and use Betti identity arguments to derive an expression for 



Acr (see ( 145 1). This is interpreted physically as the change in 
traction near the crack tip induced by the inclusion's presence; 
as such we say that the sign of Acr for any given positioning 
and configuration of the inclusion either shields or amplifies the 
propagation of the main crack. 

We conclude the paper by presenting numerical results in 
Section [9] In particular we show how <tq varies depending on 
the extent of interface imperfection and choice of material con- 
trast parameter for different loadings. These computations are 
performed for point loadings that are chosen to be illustrative of 
the suitability of our method to asymmetric self-balanced load- 
ings. We further propose a method of comparing cr with stress 
intensity factors from the analogous perfect interface problem 
and find agreement as the extent of interface imperfection tends 
towards zero. We also present computations that show the sign 
of Actq for varying location and orientation of the micro-defect. 

3. Formulation of physical and weight function problems 

3.1. Physical formulation 

We consider an infinite two-phase plane with an imperfect 
interface positioned along the positive x-axis. A semi-infinite 



Figure 1 : Geometry for the physical setup. The crack tip is placed at the origin 
of an infinite plane composed of materials with shear modulus jii occupying 
half-planes US" above and below the the crack and imperfect interface for j = 
1,2. The central point Y of a micro-defect is situated at a distance d from the 
tip of the main crack. 



crack is placed occupying the line {(x,y) : x < Q,y — 0). We 
refer to the half-planes above and below the crack and inter- 
face respectively as I1 (I) and II (2) . The material occupying Uy' 
has shear modulus [ij and mass density pj for j = 1,2. The 
anti-plane shear displacement function u satisfies the Laplace 
equation 



V 2 w(x,y) = 0. 



(4) 



The plane also contains a micro-defect whose centre is at the 
point Y; we will consider in particular elliptic inclusions al- 
though other types of defect may be incorporated into the model 
provided a suitable dipole matrix can be obtained (see for ex- 
ample Mishuris et al. (2011 1 in which micro-cracks and rigid 
line inclusions are considered). The defect g E has shear mod- 
ulus yUi n , is placed at a distance d from the crack tip, makes an 
angle (p with the imperfect interface and is oriented at an angle 
a to the horizontal as shown in Figure [T] The value of jU; n may 
be greater than or less than the value of /j. out (which may be p.\ 
or fi2 depending where the defect is placed), and so both stiff 
and soft defects can be considered. 

We assume continuity of tractions across the crack and inter- 
face, and introduce imperfect interface conditions ahead of the 
crack: 



P\ 



du 

dy 



= V2 



y=0+ 



du 
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Kfli 



du 
dy~ 



0, x > 0, 



(5) 



(6) 



y=0+ 



where the notation [m] defines the jump in displacement across 
x = 0, i.e. 

[m](jc) =«i(jc,0 + )-m 2 (jc,0"). (7) 

The parameter k > describes the extent of imperfection of the 
interface, with larger k corresponding to more imperfect inter- 
faces. We further impose prescribed tractions p± on the crack 
faces: 



Pi 



du 
dy 



P+(x), Mi 



du 
dy 



p-(x); x<0. (8) 
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These tractions are assumed to be self-balanced; that is 



u u 

J" p+(x)dx - J" p-(x)dx = 0, 



(9) 



and it is further assumed that p±(x) vanish in a neighbourhood 
of the crack tip. Although the techniques we will establish can 
be applied to any permissible loading, we will particularly fo- 
cus our attention on the case where these loadings are point 
loadings, with a loading on the upper crack face positioned at 
x - -a (where a > 0) balanced by two equal point loadings on 
the lower crack face positioned at x — -a — b and x — -a + b, 
where < b < a. This loading makes computations more dif- 
ficult to perform than for the smooth loadings considered by 



Antipov et al. ( 2001| l, but is more illustrative of the asymmetry 
of the load. 

Near the crack tip, the physical displacement behaves as 
(-iy+W AW* / (r 



TTflj 



f2 



+ (- 1V +1 (n + (- 1 ) ] e)r sin0} + 0(r A ln z r), r -> 



„2 ,„2 , 



(10) 



as demonstrated by Mishuris ( 200 1 1 > . It follows that the dis- 
placement jump is approximated by 



[kJ(jc) = k<t + 0(xln \x\), x — > 0*, 



(11) 



as the crack tip is approached along the x-axis. 

In the neighbourhood of the crack tip, the out of plane com- 
ponent of stress behaves as 

(-iy f , 

o-j lcr lnr sin 6 + cq sin 6 + (-l)Vo(w + (-l)-'0)cos G\ 

(12) 

as r — > 0, in the usual polar coordinate system and so along the 
interface, 

cr ~ cr , x — > + . (13) 

These estimates demonstrate that Fourier transforms of the dis- 
placement jump and out-of-plane stress components can be 
taken; we denote the Fourier transformation / of a function / 
by 



j f(x)e^ x dx. 



(14) 



Thus as £ — > oo, the Fourier transform of the displacement jump 
behaves as 

lum = -Ko- ir l + o(r (1+ % ^^oo, 6>o. (is) 

Moreover, along the axis, the out of plane stress component 
decays as 



& = cr^- 1 + <9(r (1+<5) ), £ -> 5 > 0. 





Unperturbed physical problem 
for imperfect interface 



Weight function configuratioi 
for imperfect interface 



Figure 2: Geometries for the unperturbed physical (Figure 2a) and weight func- 
tion (Figure 2b) setups. 



3.2. Weight function formulation 

The sought weight function U also satisfies the Laplace equa- 
tion, but with the crack occupying {(x,y) : x > Q,y — 0). We 
define the functions £y in their respective half-planes by 



dUj 

I,j{x,y) := fxj—, 7=1,2. 



(17) 



Boundary conditions analagous to the physical set-up apply. 
That is, 

Si(x,0 + ) = Z 2 (*,0"), xeR, (18) 



[f/](x)-/fE 1 (x,0 + ) = 0, x<0. 



Z!(a;,0 + ) = S 2 (x,0 _ ) = 0, x>0. 



(19) 



(20) 



We expect that along the interface, the displacement jump be- 
haves as 



[£/](*) = 0(1), x^Q-; 
lUj(x) = 0(\x\-" 2 ), x -» -co, 

while along the crack, 

[C/J(x) = c\ + C2->clog x + ct,x + o(x), x — > + , 

and 

[t/](jc) = c\ + C5 yfx + o( yfx), x — > +oo, 
where c, are constants. We further expect that 



O(l), 



0" 



0(x- 1/2 ) 



(21) 
(22) 



(23) 
(24) 

(25) 



(16) 



3.3. Derivation of Wiener-Hopf type equation for the weight 
function 

The asymptotic behaviour of Uj allows us to apply Fourier 
transforms. Moreover, the behaviour near r — demonstrates 
that the Fourier transform exists as a Cauchy value integral. Ap- 
plying the Fourier transform with respect to x 

oo 

£7/6)0 = J Uj(x,y)e^dx (26) 



and taking into account the behaviour of U at infinity, we obtain 
that the transformed solutions of Q are of the form 



Uj(t,y)=Aj(t)e-W, 



(27) 



with the corresponding expressions for tractions at y = 0* given 

by 

S^,0 ± ) = (-1)V 7 |^(^). (28) 



We define the functions by 

<D"(f) = £| v=0 +, 

0+© = [£]-*£| v=0+ . 



(29) 



(30) 



These functions are analytic in the complex half planes 

denoted by their superscripts. We expect that as £ — > oo in their 
respective domains, asymptotic estimates for ( & ± (^) are 



and near zero, 

a> + © = o© 3/2 ), 



<d-(0 = o© 1/2 ), 



0; 



(31) 



(32) 



we verify this later (see equations (55 1-(58 1). The condition of 
continuity of tractions across the crack and interface |5]) gives 
that 

mAx(g) = -/i 2 A 2 ©, (33) 

and the Fourier transform of the jump function [17] can be seen 
from d26l) to be 



=A 1 (#-A 2 (0. 



(34) 



Combining these conditions (33l-(34l, we conclude that the 



functions <t) ± (^) satisfy the functional equation of the Wiener- 
Hopf type 

(35) 



where 



<D + © = -KEfg)Q-(£), 



m = i + 



with the constant fiQ given by 



Mo 



fllfl 2 K ' 



(36) 



(37) 



This Wiener-Hopf kernel E(£) is the same as that found by An 



however different. 



4. Solution of the weight function problem 



4.1. Factorisation of the weight function Wiener-Hopf kernel 

In this section we factorise the function E(£) as defined in 
(36 1. As we just remarked, despite this function having been 
previously factorised in Antipov et al. ( 2001| l, we provide here 



an alternative factorisation which is more convenient for com- 
putations. We define an auxiliary function S» by 



.1/2,1/2 

S.(f) = tanh 



1 /2 

with the functions £ ± given by 



1 



1/2 



/=Fif. 



(38) 



(39) 



Here yf is the standard square root function with its branch cut 

1 /2 

positioned along the negative real axis. Thus £ ± are analytic 
functions in half-planes corresponding to their respective sub- 
scripts. Now, H,(£) is an even function and behaves at zero and 
infinity as follows: 



Mo 6 \/*o 



3\ 



5.(0 = 1 + |j + 0( e - 2| ^°), |f| -» oo. 



The kernel function E(£) can be factorised as 
where 



£ S,(0Ho(f), 



^1/2^.1/2 



H © = coth 



This function can itself be factorised as 
Ho© 

where 



r(i + ^) 



V 2 n/io / 



(40) 
(41) 

(42) 
(43) 
(44) 

(45) 



The functions EjJ(f) satisfy Et(£) = E (-£), with HJ© being 
regular and non-zero in the half plane Im(^) > -nfio/2. More- 
over, Stirling's formula gives that the behaviour as £ — > oo in an 
upper half plane is 

S+(0 = /3 1/2 + i/T 1/2 + y^/T 3/2 + O06- 5/2 ), f -» oo, (46) 

where j3 = i^linpo). Analogous asymptotics for S (£) are eas- 
ily obtained by noting that Eq(0 = E„~(— £). Near £ = 0, the 
asymptotics for Et© are given by 



tipov et al. (2001 1. The behaviour of the functions O*© is Ej© = 



1 21n(2)£ (tt 2 + 12 ln 2 (2))£ 2 



^ 3/ Vo 6// 2 7r 5 / 2 
The function E«© can be written in the form 

e,© = e:©e;©, f e c* 

where 



o(f), f -»0. 

(47) 
(48) 



+ 1 fin 2.(0 



(49) 



In particular, we stress that the functions E*(£) are easy to com- 
pute. Near zero, we find that 



k<& = i + — + o(f), t^o, 



where 



a 



dt, 



(50) 



(51) 



which follows from a similar derivation to that of Vellen der et| 



al. (2011 1. Moreover, behaviour near infinity in a suitable do- 
main is described by 



(52) 



These expressions again emphasise the well behaved nature 
of the functions S*(£). The 'bad' behaviour of the kernel near 
£ = is all contained in the function Eo(£) which has subse- 
quently been factorised into the product of readily computable 
analytic functions. 

4.2. Solution to the Wiener-Hopf weight function problem 
In this section we solve the Wiener-Hopf problem given in 



equation (35 i. Substituting our factorised expressions for Eg(£) 



and E*(£) into ( 35 1, we arrive at the Wiener-Hopf type equation 



1/2 



g4> + (g)£ 

ej(£)h:(£ 



-jotaioO (£>S„(£>S (£)-__ 1/2 . 



± 



(53) 



Both sides of (53 1 represent analytic functions in their respec- 



tive half-planes and do not have any poles along the real axis. 
The asymptotic estimates as £ — » oo given in ( 31 1, (46 1 and ( 52 1 



demonstrate that both sides of equation ( |53) behave as (9(1) as 
£ — > oo in their respective domains. We therefore deduce that 
both sides must be equal to a constant, which we denote M. 
We deduce that the functions ( I> ± (^) are given by 



0>+(£) 



1/2 ^O©"* ©■ 



(54) 

These expressions validate our earlier expectations (see equa- 
tions ( 3 1 1 and ( 32 1 on page |4| regarding the asymptotic esti- 
mates for <& ± . In particular, accurate estimates near zero are 
given by 



<D+(£) 



1/2 



1 (55) 



1/2 



^,0^0 + 21*2^ ^ 



7n// 



(56) 



while as ^ — » oo in the appropriate domains, 



r(i)= _^(,_!f!M to (';|, (57) 

^Kyfp^\ m £ \£/ 



£ — > oo. 



(58) 

It also follows from < [54| ) that the Fourier transform of £/ is given 
by 



/i y /OTyu E;(£)E-(^y 2 f 



1,2. 



(59) 



Expressions for the transforms of the displacement jump and 
the mean displacement across the interface are therefore respec- 
tively given by 



(60) 



(U)(g) := J (f7i(£,0 + ) + £/ 2 (f,0-)) = ^ 



2^E,-(^)E-(^y 2 ^ 



(61) 



where p„ is the dimensionless mechanical contrast parameter 

0"i -/^) 



(Mi + 



(62) 



These expressions will be useful in Section[5]where we consider 
the Betti identity in an imperfect interface setting. In particu- 
lar we note that [£/] has asymptotic expansions near zero and 
infinity as follows 



I°M>-^(i + ^*o(i 



(63) 



m. 



1/2 



(21n(2)-arW , 
1 + — — + 0(f) 

7TyU 



The function (U) behaves similarly, as 



(64) 

£ — > 00, 
(65) 



2tt 3 / 2 ^ 



1/2 



TTflQ 



(66) 



Another key difference between the imperfect and perfect in- 
terface (as considered in |Piccolroaz et al7| ( |2.009) >) cases is also 
readily seen here. Due to the condition of continuity of dis- 
placement across perfect interfaces, the function [!/](£) is a 
plus function in the perfect case, since [17] (jc) is zero for x ly- 
ing along the negative real axis. However, across an imperfect 
interface, the displacement is no longer continuous and so [{/] 
is neither a plus function nor a minus function. 



5. Betti identity in the imperfect interface setting 

In this section we refer to the physical fields for displacement 
and out-of-plane stress component as u and <x respectively, and 
to the weight function fields for displacement and stress as U 
and S respectively. We will use the reciprocal theorem (Betti 



formula) as in Willis et al. ( 1995 i to relate the physical solution 



to the weight function. 

Applying the Betti formula to the physical fields and to the 
upper and lower half plane we obtain 



CO 

/ 



{ t/(x' - x, + Mx, + ) - Hx' - x, Q + )u(x, + )) dx = 0, 



(67) 



and 



CO 

/ 



{U(x - x, Q-)cr(x, (T) - £(x' - x, 0~)u(x, - )} dx = 0. 



(68) 

These identities were proved under the assumption that the in- 
tegrand decays faster at infinity than l/R along any ray. It is 
clear from the asymptotic estimates for the physical solution 



and the weight function given in subsections 3.1 



this condition is satisfied. Subtracting ( 68 1 from ( 67 1 we obtain 



and 3.2 that 



oo 

/ 



\ U(x - x, + )cr(x, + ) - U(x - x, (T)o-(x, 0") 

[£(x' - x, + )u(x, + ) - T,(x - x, Q-)u(x, (T)]}dx = 0. 

(69) 



We split the terms for physical stress into two parts, writing 

oix, 0*) = p¥(x) + cr (+) (x), (70) 

where p± and <x (+) are defined as follows 

p ( ~\x) = ffC-JcMjc.O*), o- (+) (x) = ff(jc)cr(jc,0); (71) 

here H(x) denotes the Heaviside step function. The functions 
represent the prescribed loading on the crack faces. After 



this splitting, equation ( 69 1 becomes 

OO 

J {{Ujix' - x)o- {+ \x) - £(x' - x, 0)|wj(jc)}dx 

o 

[U{x - x, + )p + (x) - U(x - x, Q-)p-{x))dx. (72) 



7' 



We introduce notation for symmetric and skew-symmetric parts 
of the loading: 

{pt\x) = \(p ( ;\x) + pt\x)), i P f-\x) = P l ~\x) - p ( :\x). 

(73) 



This allows us to rewrite the right hand side of (72 1, giving 

oo 

J {fUjix - x)cr (+) (x) - E(x' - z,0)[w](*)}dx 



CO 

■/ 



= - {[!/](*' - x){p)(x) + (U)(x' - x)lpj(x)}dx. (74) 



We now split [17] into the sum of [t/] (±) in the spirit of (70 1, 
and similarly split [h] into the sum of [[m] (±) . We will use the 
usual notation of / * g to denote the convolution of / and g. 



Rewriting ( 74 1 using these expressions gives 



M 



(+) * rr(+) 



(-) * ^(+) . 



■ £ * hi 



l(+) . 



l(-) 



- \U\ * { P t ] - (U) * [p]<->. 



(75) 



Taking Fourier transforms in x yields 



[C/]«(0c7< + >(0 + [f/]"(^«(£) 



-2(flM w ©-2ffll"l H (fl 

= -lU}(0(p)(£) - (U}({)fpl({). (76) 

We now make use of the transmission conditions which state 
that 

[l/]<->(0 = [n]«(0 = «r«©. (77) 

This causes the second and third terms in the left hand side of 
f76) to cancel, leaving 



[£/]W©<r( + )®-X@Wl-)© 

= -[ff]©<P>(0-<tO(0W(£>. {eR. (78) 

We note that 

[P© = 1> + (^), = *"(£). (79) 

and can therefore combine the asymptotic estimates in l ] 15} , 
(fT6), d57) and d58) to yield that 



[t/](+W+>(£) 



cr i 



+ o(r (2+,5) ) 



+ 0(^ (2+,5) ) 



^ ^ oo in C + , 
(80) 

£ — > oo in C~, 
(81) 



where 5 > 0. We now multiply both sides of ( |78[ > by giving 

- (iuk&cpxo + (umym) . (82) 

Then, and similarly to the expression obtained for the perfect 
interface Betti formula approach of W illis et al.| ( |T995] >, the left 
hand side now has asymptotics at infinity (in appropriate do- 
mains) of the form 



^([C/]W(^crW(^)-S^)H(-)^)) 



1 



1 



£ + iO £ - /0 
(83) 



as £ — > oo, where the term in square brackets is the regulariza- 
tion of the Dirac delta function, namely -2m6{^). Integrating 



both sides of ( 82 1, we can arrive an expression for the constant 



cr in terms of known, readily computable functions: 



cr = 



We note that since [£/](£) and <!/>(£) behave as 0{^ 2 ) as 
£ — > oo, and the functions [/>](£) and (/?)(£) behave as bounded 
oscillations as £ — > oo for point loadings, the integrand is well 
behaved at infinity. Moreover, near £ = the integrand is also 



— 1/2 

sufficiently well behaved, acting as 0(^ + ) 



Equation ( 84 1 is a particularly important result; it gives an 



expression for the leading order of the out-of-plane component 



of stress near the crack tip (see ( 13 1) in terms of known func 



tions and acts as an imperfect interface analogue to the stress 
intensity factor from the perfect interface setting. 

5.7. The functions [p] and (p) for specific point loadings 

As we stated earlier in this paper, although the methods de- 
scribed are applicable to any permissible loading, we will later 
perform computations using the specific point loading configu- 
ration shown in Figure[T]on page|2] 

For this configuration, the loadings are defined as a point load 
on the upper crack face at x = -a balanced by two equal loads 
at x — -a + b and x = —a — b, that is 



p[ \x) - F6(x + a), p ( _ '(x) - '— (6(x + a + b) + 6(x + a - b)) . 

(85) 

The corresponding explicit expressions for (p)(x) and 
are 



F [ 1 

(p)(x) = — < S(x + a) + - (6(x + a + b) + S(x + a - b)) 



(86) 



Ipj(x) - F \5(x + a) - - (6(x + a + b) + 6(x + a - b)) \ , 



which have Fourier transforms given by 



</>>(£) = + 1)V*"** 

[p]© = -f(^-i)V i(fl+ «? 



(87) 

(88) 
(89) 



6. The unperturbed solution, u 

We will later require a method to evaluate the unperturbed 
physical solution uq and its first order partial derivatives with 
respect to x and y. This problem has been solved by Antipov 



et al. ( 2001) 1 by approximating the loading by a linear combi- 



nation of exponentials; this approximation is however not ideal 
for point loadings. 



Tractions on the upper and lower crack faces can be written 



as 



o- 1 (f,0 + )=pig) + ^g), 



(90) 



(91) 



It follows immediately from continuity of tractions across the 
imperfect interface that 

(0 = <fi(g) =: <p + (£). (92) 

We further define minus functions, <p^ and <f 2 as 

^(0 = [S]©-^ 1 (0O + ), (93) 



(94) 



We expect that the unknown functions <p + (^) and tp. (£) behave 
at infinity as 



oo, ±Im(£) > 0. (95) 



From these expressions follow the relationships 



and also 



(96) 

(97) 
(98) 



2[s](0 - 2k{o-){0 = ipi(g) + (99) 
Moreover, since transformed solutions are of the form 



ufcy) = Aj(g)e-W, 
we further have the relationships 

[S1(0=A 1 (0-A 2 (a 



and 



[a-l(0 = -|^IOuiA 1 (0+// 2 A 2 (0), 



(100) 



(101) 



(102) 



(Dig) = y(^A 2 © -piA^)). (103) 

These seven equations in eight unknowns reduce to the follow- 
ing Wiener-Hopf type equation relating ip + {^) and 

= ft(fl+|{l-^}H(f). «04) 



Noting that the term in braces on the left hand side of ( |104| > is 
the function we earlier defined as E(£) and have aheady suitably 
factorised, we can write 

-*H(^ + (£) - KE(f)(p)(.f) = <p[(g) + KA(£)$pl(g), (105) 



where 



A(0 = I/l_^ 
? 2 \ |f| 



Recall that H(£) can be factorised in the form 
B(0 = nfi B + (Z)B-(£), 



(106) 



(107) 



where we have defined the functions B ± (^ r ) for the sake of nota- 
tional brevity by 



1/2 



fl-(g)= B ° ( ^' (g) , .108, 



which are analytic in the half planes indicated by their super- 
scripts. These functions have behaviour near zero and infinity 
given by 

B*(j~) = 0(g~ 1/2 ), £^0, (109) 



Thus 



= 0(1), 



% — > oo. 



(111) 

We can decompose the final term on the right hand side as usual 
into 

o (?) 

where L ± (^) are given by 



DO 



S-08)' 



(113) 

for ^ e C*. We expect that behave as O^ 1 ) as £ -> oo. 
The Wiener-Hopf equation becomes 



o (?) 



(114) 



Both terms on each side of (114i decay as 0(1/%), £, — > oo. 
Moreover, each side is analytic in the half-plane denoted by 
the superscripts. Liouville's theorem yields that both sides are 
equal to zero, and so 



<p\(?) = L-(g)BT(g). 



(115) 



These expressions verify that our expectations of the behaviour 



of tp + {%) and <p 1 (g) as ^ — > oo were correct. Moreover, (98i 
enables us to express (p^S) as 



(116) 



Condition (99 1 then yields an expression for the transform of 



the displacement jump 

M/ » <pj(& 

[«](£) = - J 2- + -\- + K(<r)(g) 

= + ^> + (£ + + K - [p](0, (117) 



from which we can obtain expressions for Aj(£) and AaCf) as 
follows 



Am = — ^ {? + «o + <pm + \ 



(118) 



These expressions now enable us (see ( | 1 00[ >) to compute the 
Fourier transform of the unperturbed solution (i.e. the setup 
with no small defect present) Uj(£, y) for any y. 

7. Perturbation analysis 

We shall construct an asymptotic solution of the problem us- 



ing the method of |Movchan et al. ( 1995[ ), that is the asymptotics 
of the solution will be taken in the form 



(1 10) u u2 (x, e) = ufl(x) + eW (1) (£) + e 2 u ( -\(x) + o(s 2 ), 



In (120 1, the leading term ufl(x) corresponds to the unper- 
turbed solution, which is described in the previous section. The 
term eW (1) (£) corresponds to the boundary layer concentrated 
near the defect and needed to satisfy the transmission condi- 
tions for the elastic inclusion g B 

du m du ou i 

Win = «out, f*in.-TT~ = f*oat-Z— on d g s . (121) 

on on 

The term s u\ Jx) is introduced to fulfil the original bound- 
ary conditions (4) on the crack faces and the interface condi- 
tions (2), (3) disturbed by the boundary layer; this term, in 
turn, will produce perturbations of the crack tip fields and cor- 
respondingly of the constant cr () . 

We shall consider an elastic inclusion, situated in the upper 
(or lower) half-plane. The leading term ufl clearly does not 



e -> 0. 
(120) 



satisfy the transmission conditions ( 121 1 on the boundary dg s . 
Thus, we shall correct the solution by constructing the boundary 
layer W (1) (£), where the new scaled variable £ is defined by 



x-Y 



(122) 



with Y = (X, Y) being the "centre" of the inclusion g s (see Fig- 
ure [TJ. 

For W (I) (£) = {WP, £ e g; W^ v £ e R 2 \ g} we consider 
the following problem 

VX 1 ^) = 0, €*g, V 2 W^ t (£) = 0, f e R 2 \ g, (123) 



where 



= s- 1 g e = {^€R z :s€ + Yeg £ }. 



The function remains continuous across the interface dg, 
that is, 



W 



(i) 



and satisfies on dg the following transmission condition 



dn 



dn 



(124) 



as e — > 0, where n — is an outward unit normal on dg. The 
formulation is completed by setting the following condition at 
infinity 



W. 



(i) 



as |f | — > oo. 



(125) 



The problem above has been solved by various techniques and 



the solution can be found, for example, in Movchan et al. 
([1995). 

Since we assume that the inclusion is at a finite distance from 
the interface between the half-planes, we shall only need the 
leading term of the asymptotics of the solution at infinity. This 
term reads as follows 



M w 



+ 0(|fr 2 ) as Z^co, 



where Ai is a 2 x 2 matrix which depends on the characteristic 
size ( of the domain g and the ratio fi out //J.i n ; it is called the 
dipole matrix. For example, in the case of an elliptic inclusion 
with the semi-axes t a and £/, making an angle a with the positive 
direction of the x-axis and y-axis, respectively, the matrix At 
takes the form 



where 



M 



1 + cos 2a 1 
+ 



■tt(i + *-i)s, 



(127) 



e + //* 1 + e/j* 

(1 - e)(/K* - 1) sin2o' 



■cos 2a (1 - e)(//* - 1) sin 2c 

(e+/**)(l +<?//*) 
1 - cos 2a 1 + cos 2a 



(e + yu*)(l + e//*) e + /j* 1 + e/j* 

(128) 

e = tblla and fj.+ = jU ou t/A<in- We note that for a soft inclusion, 
y«out > /"in, the dipole matrix is negative definite, whereas for a 
stiff inclusion, /i out < /jj n , the dipole matrix is positive definite. 
In the limit /i m — > oo, we obtain the dipole matrix for a rigid 
movable inclusion. In the case of an elliptic rigid inclusion, we 
have 

(129) 



where 



Atrfg = + l)S rig , 

h+(a) + eh-(a) (l-e)sin2o' 
(1— e)sin2a h-{a) + eh+{a) 



Here we have defined the functions h ± (a) — 1 ± cos 2a for 
brevity of notation. The term sW (l) (t;) in a neighbourhood of 
the x-axis written in the x coordinates takes the form 



eW (1) (f) = e 2 w m (x) + o(e 2 ), s 0, (131) 



where 



w ( V) = -^[v,«H y ] 



M 



x-Y 



(132) 



As a result, one can compute the average e 2 (cr (1) ) and the 
jump £ 2 Jcr (I) ] of the "effective" tractions on the crack faces in- 
duced by the elastic inclusion g s . Since du. Jdy = -<9w (1) /dy 



must hold on the crack line (to satisfy the original boundary 
conditions (4)), this gives for x < 



(cr a) ) (-) (jc) = _i (Ml+iU2 )^ :=j p(-) W) 

2 dy 



where 



oy 



(133) 
(134) 



dy 



i* - in 2 

77- L ITJ r-YH 



(135) 



Additionally, we can compute the transmission conditions for 



(126) 

the functions i/^ across the interface. In order for the perturbed 



solution U\2 in * 120 1 to satisfy the original transmission condi- 



tions (2) and (3), the following relations must hold for x > 



(i>H(+V 



(x) = ^ (1) ) (+) (x)+^0u 1+yU2 ) 



2 ri n ' dy 
: = k<o" (1) > (+) (x) - kP (+) (x), 



3y 



(136) 
(137) 



8. Model problem for the first order perturbation 

The constant ctq which describes the traction near the crack 
tip (see ( fT~3] >) is expanded in the form 

o-o = of + e 2 Acr + o(s 2 ), e -> 0. (138) 

Our objective is to find the first order variation A<x () . 

Let us consider the model problem for the first order per- 
turbation u w and write the corresponding Betti identity in the 
form 



f 



1 ( 

[£/](*' - x)(cr w )(x) + (U)(x' - x)[cr (1) ](x) 

-<E>(x'-x)[ M (1) ](x)ldx = 0. 



(139) 



(130) This follows immediately from (69 1 by noting that [S] = 0. We 



split the terms for stress into two parts, 



Am 



Q (+) + Q' 



(140) 



observing that in contrast to the zero order problem where the 
load is described by (73 i, the terms with superscript (+) are non- 
zero since the presence of inclusions induces stresses along the 
imperfect interface and should be taken into account. Equation 
(jT39j becomes 

[lU}(x' - x)< ( r (1) > (+) (x) - <E>(x' - x)[ M (1) ](x)} dx 

OO 

= - f Wl& - x)P { -\x) + (U)(x' - x)Q ( -\x) 

+(U)(x' -x)g (+) (x)ldx. 



We now split [17] into the sum of [t/]* and similarly split [«] 
into the sum of [«]*. This gives 

[t/] (+) * <cr (1) ) (+) + [C/JH * <<r«> (+ > - <Z> * [m (1) ] (+) 
-<2> * [w (1) ] H = -Wl * J* - <t/> * 6 M - (C/> * e (+) . 



(141) 



Taking the Fourier transform in x yields 

-<I>(^)Im (1) ] + (^) - <2)(^)[i7 (1) ]-(^) 

= -$Uj(g)T(g) - (U){m^) - (U)({jQ + ({). 

We now make use of the transmission conditions 



(142) 



[T7]-(# = *<£>(£), 



7(1)H + / 



© = ^ 1) > + ©-^'(a 



(143) 



thus obtaining 



(144) 



The same reasoning used in Section [5] allows us to derive the 
integral representation for A<x in the form 



Ao-n 



[^<2>-(^)p + (^) + ^<i/)(^)e + (^] df 

DO 



(145) 

This important constant has an immediate physical meaning. If 
Actq = then the defect configuration is neutral; its presence 
causes zero perturbation to the leading order of tractions at the 
crack tip. Otherwise, if Acr < 0, the presence of the defect 
causes a reduction in the crack tip traction and so shields the 
crack from propagating further. Finally, if A<x > then the 
defect causes an amplification effect and so can be considered 
to be encouraging the propagation of the main crack. 

9. Numerical results 

9.1. Computations of o~o 

In this section we present results of computations obtained 
by following the methods previously described in this paper. 
All results have been computed using MATLAB. 

Figure [3]plots cr against p.*, showing how the constant from 
the asymptotic expansion at the crack tip cr varies with differ- 
ently contrasting stiffnesses of materials. Recalling that 



(146) 



we note that when p* is near to -1, this corresponds to p2 s> 
p\. That is, the material occupying the region below the crack 




-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



Figure 3: Plot of <tq against Both cases plotted here use the parameters 
k, = 1 and a = 1 , but with different values for b, which controls the separation 
between the point loadings. The red plot has b = 3/4 while the blue plot uses 
6=1/4. 



is far stiffer than the material above the crack. As this limit 
is approached, the precise locations of the point loadings on 
the lower face of the crack decrease in importance, since the 
material becomes sufficiently stiff for the material to act as an 
almost rigid body; this explains the meeting of the two lines at 
H» = -1. 

In Figure [4] we present a log-log plot of <r against k„, the 
dimensionsless parameter of interface imperfection defined as 
k* = Ar(yUi + ii-i)la. This has been computed for different values 
of //* (describing the contrast in material stiffnesses) and also 
for different values of b (describing the separation distance be- 
tween the point loadings) while keeping a fixed (a = 1). The 
solid lines correspond to b = | while dotted lines represent 
b = 4 and different colours correspond to different values of 
fi t : green corresponds to //„ = -0.8, blue to fj.* = and red to 
[£» = +0.8. 

Bearing in mind our remarks regarding Figure [3] we would 
expect that changing the value of b would have the greatest im- 
pact for values of //„ near +1. This is indeed the case in Figure 

El 

Also plotted in Figure [4] is a grey dotted line that is tangent 
to the curves (which run parallel) as K t — > 0; this tangent has 
slope -i indicating that <x = 0(k~ l/2 ) as k* — > 0. As k* — > 0, 
the interface becomes almost perfect, and so the square-root 
behaviour associated with fields near crack tips in the perfect 
interface setting is not unexpected. Moreover, as k — > +oo, the 
curves on the log-log plot have slope -1, implying that cr = 
0(k~ 1 ) as K t — > +oo. 

9.2. Comparison of o~q with stress intensity factors from the 
perfect interface case 

In this subsection we discuss an approach which enables a 
comparison to be made between imperfect and perfect interface 
situations. 

Comparing the fields directly is not a simple task since in 
the perfect interface case the stresses become unbounded at the 
crack tip, exhibiting asymptotic behaviour of <x = 0(r~ l t 2 ), 
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Figure 4: Log-log plot of <xo against k, for differently contrasting materials. 



Figure 6: Plot of the sign of Ao"o for varying a and <p with the problem param- 
eters as described in Figure. The darker shaded areas are those (0, a) for which 
Aero > while paler regions have Actq < 0. 



1 

0.99 
0.98- 
0.97 
&. 0.96 
0.95 
0.94 
0.93 
0.92 





-W a =aa 
-(«.)>= as 

-W> = o.8 



0.2 



Figure 5: Plot of the ratio r as defined in (149) for four different values of Qi,)2- 

r — > 0. In the imperfect setting, we have derived the leading 
order of stresses at the crack tip, cr , which is independent of 
r. Moreover, different normalisations may make comparisons 
difficult. 

However, given two particular pairs of materials with con- 
strast parameters (/j*)i and (/i»)2 say, we might expect the 
dimensionless ratios of stress intensity factors (^i])i/(.^jnh 
(from the perfect interface case) and (cro) 1/(1x0)2 (imperfect 
case) to be similar for small k*. In the perfect interface case, 
the stress intensity factor (derived in Piccolroaz et al. ( 2012| l) is 
given by 

*$ = -VI jV> ( - r > + Y m - r) } rU2dn (147) 



As derived earlier in Section [5] the leading order of tractions 
near the crack tip in the imperfect interface case is given by 



Figure [5]plots the ratio 



r(K t ) 



(o-o(^)) 1 /(o-q(a- ) ,))2 
( K iun/( K iiih 



(148) 
(149) 



for < k* < 1 with (yUt)i =0 fixed and for four different values 
of (yi,)2. The loadings used are balanced; a point loading on the 
upper crack face at x = -1 is balanced by two equal loadings at 
x = -1.25 and x = -0.75. 

We see from the plot that as k* — > 0, r(K*) — > 1. This pro- 
vides some verification of the accuracy of our computations and 
demonstrates that the comparison of ratios approach for small k 
again the perfect interface case is useful. 

9.3. Computation of Actq 

We now present numerical results for the perturbed problem 
computed using MATLAB. Figure|6]shows the sign of Aero for 
a specific configuration. To reduce the computational task here, 
we have used smooth loadings with the tractions on the upper 
and lower crack faces of the form 



P+(x) 



P-(x) 



x<0. 



(150) 



and the imperfect interface has k* = 1. The inclusion is stiff, 
with the contrast between the internal and external materials of 
the inclusion given by //* =5. 

The figure clearly shows the regions for which crack growth 
is encouraged or discouraged for this configuration. How- 
ever, we make the observation that different analysis should 
be sought when <p is particularly close to zero since this cor- 
responds to the crack being placed near the imperfect interface 



which contradicts the assumption made before equation ( 126 1. 
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